Large-scale computer simulation of fully developed turbulent channel flow with heat transfer

Recently, with the advent of supercomputers, there has been considerable interest in the use of direct numerical simulation to obtain information about turbulent shear flow at low Reynolds number. This paper presents a pseudospectral technique to solve the full three-dimensional time-dependent Navier-Stokes and advection-diffusion equations without the use of subgrid-scale modelling. The technique has not been previously used for fully developed turbulent channel flow simulation and is based on methods applied in other contexts. The emphasis of this paper is to provide a reasonably detailed account of how the simulation is done rather than to present new calculations of turbulence. The details of an algorithm for turbulent channel flow simulation and the grid and time step sizes needed to integrate through transient behaviour to steady state turbulence have not been published before and are presented here. Results from a Cray-2 simulation of fully developed turbulent flow in a channel with heat transfer are presented along with a critical comparison between experiment and computation. The first- and second-order moments agree well with experimental measurements; the agreement is poor for higher-order moments such as the skewness and flatness near the walls of the channel. Detailed information given about the effects of spatial grid resolution on a computed results is important for estimating the size of the computation required to study various aspects of a turbulent flow.     

Pulsed-Field-Gradient Measurements of Time-Dependent Gas Diffusion

Pulsed-field-gradient NMR techniques are demonstrated for measurements of time-dependent gas diffusion. The standard PGSE technique and variants, applied to a free gas mixture of thermally polarized xenon and O₂, are found to provide a reproducible measure of the xenon diffusion coefficient (5.71 × 10⁻⁶m²s⁻¹for 1 atm of pure xenon), in excellent agreement with previous, non-NMR measurements. The utility of pulsed-field-gradient NMR techniques is demonstrated by the first measurement of time-dependent (i.e., restricted) gas diffusion inside a porous medium (a random pack of glass beads), with results that agree well with theory. Two modified NMR pulse sequences derived from the PGSE technique (named the Pulsed Gradient Echo, or PGE, and the Pulsed Gradient Multiple Spin Echo, or PGMSE) are also applied to measurements of time dependent diffusion of laser polarized xenon gas, with results in good agreement with previous measurements on thermally polarized gas. The PGMSE technique is found to be superior to the PGE method, and to standard PGSE techniques and variants, for efficiently measuring laser polarized noble gas diffusion over a wide range of diffusion times.     

On a new system of nonlinear A-monotone multivalued variational inclusions

In this paper, we introduce and study a new system of nonlinear A-monotone multivalued variational inclusions in Hilbert spaces. By using the concept and properties of A-monotone mappings, and the resolvent operator technique associated with A-monotone mappings due to Verma, we construct a new iterative algorithm for solving this system of nonlinear multivalued variational inclusions associated with A-monotone mappings in Hilbert spaces. We also prove the existence of solutions for the nonlinear multivalued variational inclusions and the convergence of iterative sequences generated by the algorithm. Our results improve and generalize many known corresponding results.     

Development of a finite-element analysis methodology for the propagation of delaminations in composite structures

Analysing the collapse of skin-stiffened structures requires capturing the critical phenomenon of skin-stiffener separation, which can be considered analogous to interlaminar cracking. This paper presents the development of a numerical approach for simulating the propagation of interlaminar cracks in composite structures. A degradation methodology was introduced in MSC.Marc, which involved the modelling of a structure with shell layers connected by user-defined multiple-point constraints (MPCs). User subroutines were written that employ the virtual crack closure technique (VCCT) to determine the onset of crack growth and modify the properties of the user-defined MPCs to simulate crack propagation. Methodologies for the release of failing MPCs are presented and are discussed with reference to the VCCT assumption of self-similar crack growth. The numerical results obtained by using the release methodologies are then compared with experimental data for a double-cantilever beam specimen. Based on this comparison, recommendations for the future development of the degradation model are made, especially with reference to developing an approach for the collapse analysis of fuselage-representative structures. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 43, No. 1, pp. 15–42, January–February, 2007.     

Spectral shift function and trace formula

The complete proofs of Krein’s theorem on the spectral shift function and the trace formula are given for a pair of self-adjoint operators such that either (i) their difference is trace-class or (ii) the difference of their resolvents is trace-class. The proofs, essentially due to Krein, is based on Herglotz’s theorem on the boundary value of the analytic functions whose imaginary part is non-negative on the upper half plane, and an almost optimal class of functions are obtained for which the trace formula is valid. Also an alternative method based on Weyl-von Neumann’s theorem for self-adjoint operators, avoiding the complex function theory and inspired by Voiculescu’s work, is given for the first case. Furthermore, some applications of the spectral shift function have been discussed.     

聚合物纳米孔隙增透膜制备工艺的研究

论述了聚合物纳米孔隙增透膜的制备工艺流程,分析了聚合物材料分子量、实验环境温度和湿度、溶剂挥发性等条件对纳米孔隙增透膜的影响。研究表明,聚合物材料分子量的增大、温度的降低、湿度的升高以及采用挥发性弱的溶剂都将导致增透膜孔隙尺寸的增大,孔隙越大其对光的散射损耗就会增大,所以增透膜的透过率就越低。通过大量的试验分析得出一组较理想的工艺参量:使用低分子量的聚合物材料(小于15 kg/mol),环境温度大于25℃、环境相对湿度小于30%,在采用低沸点的溶剂如四氢呋喃等措施下可有效降低增透膜散射损耗。     

An improved linearization strategy for zero-one quadratic programming problems

We present a new linearized model for the zero-one quadratic programming problem, whose size is linear in terms of the number of variables in the original nonlinear problem. Our derivation yields three alternative reformulations, each varying in model size and tightness. We show that our models are at least as tight as the one recently proposed in [7], and examine the theoretical relationship of our models to a standard linearization of the zero-one quadratic programming problem. Finally, we demonstrate the efficacy of solving each of these models on a set of randomly generated test instances.     

Comparison between Eulerian and Lagrangian semi-analytical models to simulate the pollutant dispersion in the PBL

In this work we display a numerical comparison, under statistical and computational point of view, between semi-analytical Eulerian and Lagrangian dispersion models to simulate the ground-level concentration values of a passive pollutant released from a low height source. The Eulerian approach is based on the solution of the advection–diffusion equation by the Laplace transform technique. The Lagrangian approach is based on solution of the Langevin equation through the Picard’s Iterative Method. Turbulence inputs are calculated according to a parameterization capable of generating continuous values in all stability conditions and in all heights of the Planetary Boundary Layer (PBL). Numerical simulations and comparisons show a good agreement between predicted and observed concentrations values. The comparison reveals the main advantages and disadvantages between the models.     

On the evolution of two-dimensional patterns in an inviscid liquid sheet

We perform an analysis of the pattern formation for a moving sheet of inviscid fluid. The sheet, which is assumed to have an infinite horizontal extent, moves at some prescribed velocity into a passive surrounding gas. The sheet’s thickness is assumed much smaller than the horizontal scale of the fluid motion. By considering a system that is symmetric with respect to the horizontal planes, long scale asymptotics are used to reduce the full governing equations in three dimensions to a set of three coupled nonlinear partial differential equations for the horizontal components of the velocity field and the height of the interface profile. The interfacial conditions consisting of the kinematic and normal stress balance are incorporated into these evolution equations. Investigations are carried out as function of the sole dimensionless parameter, namely the Weber number. A small amplitude stability analysis around the planar gas–liquid interface reveals that wave patterns in the form of traveling plane waves occur subcritically, and are therefore unstable. The reduced evolution equations are solved numerically for fixed values of the Weber number. Since the reduced system of equations is homogeneous, the wave motion is generated by initial conditions. Five initial conditions have been imposed: one-dimensional rolls, two-dimensional squares, two-dimensional hexagons, two-dimensional ridges, and smooth peaks. The ensuing evolution of the liquid sheet’s shape and corresponding flow fields are described by illustrations of the changes in the sheet’s morphology with time.     

13.5—14.6MeV中子能区锗的同位素反应截面的测量

报道了在13.5—14.6MeV中子能区用活化法测得的⁷⁶Ge(n,2n)⁷⁵Ge, ⁷⁰Ge(n,2n)⁶⁹Ge, ⁷⁰Ge(n, p)⁷⁰Ga, ⁷²Ge(n,p)⁷²Ga, ⁷³Ge(n,p)⁷³Ga, ⁷²Ge(n,α)^69mZn和⁷⁴Ge(n, α)^71mZn的反应截面值. 中子注量用⁹³Nb(n,2n)^92mNb反应截面得到. 单能中子由T(d,n)⁴He反应获得. 同时还列举了已收集到的文献值以作比较.